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Dynamical Triangulation with Fluctuating Topology 
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We consider a dynamical triangulation model of euclidean quantum gravity where the topology is not fixed. 
This model is equivalent to a tensor generalization of the matrix model of two dimensional quantum gravity. A 
set of moves is given that allows Monte Carlo simulation of this model. Some preliminary results are presented 
for the case of four dimensions. 



1. Introduction 

In the dynamical triangulation model of quan- 
tum gravity one usually considers systems where 
the spacetime topology is fixed (although the 
topology of space is allowed to change with time). 
In a path integral formulation of quantum grav- 
ity, summing over the topologies seems a natural 
thing to do. 

Because the typical curvature fluctuations be- 
come larger at smaller scale, allowing an arbitrary 
topology will result in a very complicated struc- 
ture at the Planck scale [ |l|. Such a spacetime 
which is full of holes is commonly called space- 
time foam. It was estimated by Hawking [ g 
that the dominant contribution to the path inte- 
gral would come from spacctimes where the Eulcr 
characteristic \ is of the order of the volume of 
the spacetime in Planck units. 

2. Definition of the model 

The partition function of the model in d dimen- 
sions is 



Z(n,N d )= J2 cxp0dV d _2). 

T(N d ) 



(1) 



This expression is the same as for a fixed topology, 
but here the sum is over all possible ways to glue 
a fixed number Nd of d-simplices together, main- 
taining orientation (i.e. only identify d— 1 dimen- 
sional faces with opposite orientation). Because 
the faces have to be glued in pairs, the number of 
simplices Nd must be even in odd dimensions. 



At this stage we include disconnected config- 
urations. Because the action is a sum of the 
actions of the connected components, the Boltz- 
mann weight factorizes, which means that the lo- 
cal physics will not change. However, including 
these configurations will tell us something about 
the chance to obtain a particular size of connected 
component. A universe which is most likely to be 
split up into many small parts seems an unlikely 
candidate for the real world. 

Naively gluing together simplices will not re- 
sult in manifolds, but only in pseudomanifolds. 
Several ways to deal with this problem are con- 
ceivable. First, as nobody knows what spacetime 
looks like at the Planck scale, one could argue 
that this is not a problem. Second, one can (at 
least for d < 4) locally deform the resulting pseu- 
domanifold to turn it back into a manifold by re- 
moving a small region around the singular points 
and pasting in a regular region. This can be done 
while changing the total curvature only by an ar- 
bitrarily small amount. 

Third, these configurations might be unimpor- 
tant in the limit k — > oo (which will have to 
be taken, see below). E.g. in three dimensions 
for each fixed Nq (and volume N3) the number 
of edges (which couples to n) is maximal if and 
only if the configuration is a real (i.e. non-pseudo) 
manifold. See [ [I| for an explanation. This is less 
clear in four dimensions, though, because a sim- 
ilar reasoning only results in a closed neighbour- 
hood of a point being bounded by a simplicial 
manifold (as opposed to a pseudomanifold) , but 
not necessarily a sphere. 



At low k the connected configurations will con- 
tribute most as they have the highest entropy. At 
high k the disconnected configuration will con- 
tribute most as they have the lowest action. It is 
a priori not clear whether this change occurs grad- 
ually or whether there is a phase transition. Sup- 
pose for a moment that there is a sharp crossover 
at some k c depending on the volume. Because the 
number of connected configurations rises faster 
with the volume than the number of disconnected 
ones, the value of k c will increase logarithmically 
with the number of simplices. This means that 
in a possible continuum limit, the value of k will 
have to be taken to infinity. 

3. Tensor model 

We can formally write down a tensor model 
which is a generalization of the well known one 
matrix model of two dimensional quantum grav- 
ity (see e.g. [|| and references therein). The par- 
tition function of this tensor model is written in 
terms of a fc-dimcnsional tensor M of rank d. The 
tensor M is invariant under an even permutation 
of its indices and goes to its complex conjugate 
under an odd permutation of its indices. The par- 
tition function for three dimensions is 



9Mabc M ade M bdf M cef j ■ 



and for four dimensions it is 



(2) 



(3) 
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Using the properties of M , the action can easily 
be seen to be real. The generalization to more 
dimensions should be obvious from these expres- 
sions. 

These expressions are only formal, because the 
interaction term can be negative for any g and is 
of higher order than the gaussian term. There- 
fore, these integrals will not converge. However, 
if we expand these expressions in g, each term of 
the expansion is well defined and the dual of each 



of its Feynman diagrams is a DT configuration of 
size N d equal to the order of g used. Each prop- 
agator carries two sets of d indices. Each vertex 
of the tensor model corresponds to a d-simplex, 
and each propagator between them corresponds 
to the identification of two d— 1 dimensional faces 
of simplices. The sets of indices at each end of a 
propagator must be an odd permutation of each 
other, making sure that the simplicial complex 
has an orientation. The contribution of a partic- 
ular diagram is 



gir* k if*-* = exp(\n(g)N d + \n(k)N d _ 2 ), 



(6) 



which is precisely proportional to the Boltzmann 
weight of the DT configuration according to the 
Regge-Einstein action, with the identification k = 
Ink. 

A similar model using hypercubes was intro- 
duced in [||]. The model described above has been 
discussed for three dimensions in [ [|] . A different 
generalization of the matrix model where the di- 
mension of the matrix couples to the number of 
points in the DT configuration (which means one 
does not get the Regge-Einstein action in more 
that 3 dimensions) has been discussed in [ o . 

4. Monte Carlo simulation 

A move can most easily be described in the ten- 
sor model formulation. It consists of first cutting 
two propagators and then randomly reconnect- 
ing them. Due to the orientability, there are dl/2 
ways to connect two propagators. In the dynam- 
ical triangulation model, this corresponds to cut- 
ting apart the simplices at each side of two of the 
d — 1 dimensional faces and pasting these faces 
together. 

Unlike the case of fixed topology with (k,l) 
moves, one can easily see that these moves are er- 
godic. The number of moves needed to get from 
one configuration to another is 0(N d ), raising 
none of the computability problems associated 
with the fixed topology [ fij. The non-existance 
of a classification of four-topologies and their un- 
recognizability is usually mentioned as a problem 
for the summation over topologies. In dynami- 
cal triangulation, however, it seems to be more a 
problem for fixing the topology. 
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Figure 1. The number of connnected components 
as a function of k in the four dimensional model. 



We use the standard Metropolis test to accept 
or reject the moves. Because (again unlike the 
fixed topology case) the number of possible moves 
does not depend on the configuration, detailed 
balance is easily obtained. 

One could restrict the simulation to connected 
configurations by checking connectedness for each 
move accepted by the Metropolis test. This might 
be rather slow, because this is not a local test. 
Also, although it seems very plausible, it is not 
clear whether this would be ergodic in the space 
of connected configurations. 

5. Results 

We have simulated this model in four dimen- 
sions at a volume JV4 = 500. We used a hot start, 
that is a configuration with completely random 
connections. This would be an equilibrium con- 
figuration at k = 0. 

The number of connected components is plot- 
ted in figure |[ We see that at low k the av- 
erage number of components is almost one, while 
at high k the average number of components is al- 
most equal to the maximum number of Nd/2. Al- 
ready at this low number of simplices, the change 
between few and many components looks quite 
sharp. 

Unfortunately, the acceptance rates for these 



moves are quite low for k near k c , of the order of 
0.1%. One of the reasons for this is that the pro- 
posed moves are not local in the sense that they 
can connect any two points in the complex. This 
means that simulating much larger systems will 
probably not be feasible with only these moves. 

6. Discussion 

The main problem to investigate is the exis- 
tence of a sensible continuum limit. Although, 
due to the factorially increasing number of con- 
figurations, the grand canonical partition func- 
tion in the normal definition does not exist, this 
does not exclude that the local behaviour of the 
system might show scaling for large volumes. 

In two dimensions a limit of the grandcanoni- 
cal partition function, the so-called double scaling 
limit, is known [ffl. In this limit we also see that 
the coupling k has to be taken to infinity. 
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